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Abstract 

A two species reaction-diffusion modei, in which particies diffuse on a 
one-dimensionai iattice and annihilate when meeting each other, has been 
investigated. Mean field equations for general choice of reaction rates 
have been solved exactly. Cluster mean field approximation of the model 
is also studied. It is shown that, the general form of large time behavior 
of one- and two-point functions of the number operators, are determined 
by the diffusion rates of the two type of species, and is independent of 
annihilation rates. 
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1 Introduction 



Recently properties of systems far from equilibrium have been studied by many 
people. Although mean field techniques may give qualitatively correct results 
for higher dimensions, for low-dimensional systems fluctuations have impor- 
tant roles. Different methods have been used to study reaction-diffusion mod- 
els, including analytical and approximational methods. Among them, there 
are methods to obtain some quantities can be obtained exactly. For example 
in mUCIl, imposing some constraints on the reaction rates leads to a closed 
set of equations for average number densities in each site. The empty inter- 
val method, is another method, which has been also used to analyze the one 
dimensional dynamics of diffusion- limited coalescence • The most gen- 

eral one-dimensional reaction-diffusion model with nearest-neighbor interactions 
that can be solved exactly through empty interval method, has been introduced 
in [Sj. Empty interval method has been also generalized in Different 
methods has been introduced to calculate different quantities exactly. However, 
exactly solvable models are only a special class of reaction-diffusion models, and 
so people are motivated to use also approximate methods to understand the role 
played by fluctuations. In ^I] a two species model has been considered. In this 
model there are three competing reactions A A — > 0, BB — 0, and AB — 0. 
Asymptotic density decay rates of the two type of species for a special choice of 
parameters have been studied using the Smoluchowski approximation and also 
field theoretic renormalization group techniques. A similar model focusing on 
the same diffusion rates for the two type of species has been studied in ^21- Field 
theoretic renormalization group analysis suggest that contrary to the ordinary 
mean-field technique, the large time density of the minority species decays at 
the same rate as the majority ones in one-dimensional case. Although ordinary 
mean-field technique, generally do not give correct results for low-dimensional 
systems, its generalizations such as cluster mean-field may give correct results. 
Any how, in the mean field approximation at most one-point functions may be 
obtained. To obtain more-point functions one should use other methods. One 
possible way is, to use a generalization of mean field known as the cluster mean 
field approximation. 

One of the topics, which have absorbed many interests in recent years, is 
non-equilibrium phase transitions. There are examples, in which mean field 
(MF) solutions are not valid, but its generalization, cluster mean field (CMF) 
gives qualitatively correct results (13114115116) . A coagulation-production model 
is recently considered in Although MF equations do not give correct re- 

sults, CMF approximation predicts phase transition, supported also by Monte 
Carlo simulations. Steady state properties in the absorbing phase of Id pair 
contact process model are also investigated using Monte Carlo simulations and 
the cluster approximation. The cluster approximation qualitatively confirms 
the numerical results |14j . 

The scheme of the paper is as follows. In section 2, The mean field equations 
for general parameters have been solved exactly. It is seen that, the large 
time behavior of the average densities depend both on initial average densities 



f 



and reaction rates, and are independent of diffusion rates. In section 3, The 
cluster mean field equations for one- and two-point functions have been solved 
numerically. It is shown that the general large time behavior is determined by 
the diffusion rates. 



2 The Mean Field Approximation 

The model addressed in this article is a two-species exclusion reaction-diffusion 
model. That is, each site is a vacancy (0) or at most occupied by a particle A 
or B. The interaction is between nearest sites, and the adjacent sites interact 
according to the following interactions with the indicated rates. 

A9 <-» %A D A 

B% <-» 0B D B 

AA -> 00 A/2 

BB -> 00 A'/2 

AB -> 00 8/2 

BA — >• 00 5/2, (1) 

We consider translationally invariant initial conditions. In the mean-field ap- 
proximation, diffusion rates do not have any effect on the evolution equations 
of average number densities. The mean-field equations for the average densities 
a := (A) t and b := (B) t are 

— = -Aa 2 - Sab 
dt 

^ = -X'b 2 - Sab. (2) 
dt y ' 

The large time behaviors of these equations for special choices of parameters 
have been studied in |12lll| . Now, we want to solve these equations exactly and 
then we will show that there are cases which are not considered in |12llllj , and 
give qualitatively correct result for large time behaviors, although the exponent 
of the decay rate is not correct. 
Consider the following cases. 



I) A = A'. 

The evolution equation for u := b/a, is 



^ = (X-5)u(l-u)a. (3) 



Using (|2I3|) . it is seen that 

du (A — S)u(l — u) 



da (X + Su)a ' ^ 
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which can be integrated to 



u - 1 



1+8/X / \ 1-S/X 

(^) = (-) , (5) 
u \a J 



where uq and ao are the initial values of u and a, respectively. Now we can 
obtain the large time behavior of the average densities. It is seen that the large 
time behavior of u depends on the ratio 5/ A. 

LI) 6 > A 

At large times, obviously a — ► 0, so it is seen from (0 that depending on the 
initial value uo, two case may occur 

At large times u — * oo =>- u ~ a A / 5_1 , ~ a A / <5 
At large times u — > => u^a 5 ^^ 1 , b~a s l x 

Assuming an imbalance in the initial average densities, for example ao > bo 
(uq < 1), (J2J gives the large time behavior of a(t), and u(t) as 

a(t) - r 1 

u^-t 1 - 5 ^, (6) 

which means that for S > A, in the mean-field approximation the minority 
species dies out earlier than the majority one, and the decay exponent of u(t) 
is independent of diffusion rates. 

1.2) 6 < A 

As a consequence of the large time behavior of a, a — > 0, it is seen from (JSJ, 
that at large times u — ► 1. Defining e :=\ 1 — u |, 

l-S/X 

e ^ a T+m. (7) 

To obtain the large time behavior of a and u, we should use again J^IIHJ^ which 
give 

a(t) - r 1 

| 1 - u(t) |~ t - Trffx ; (8) 

which means that at large times both the minority and the majority species 
decay with the same rate. The exponent of decay rate does not depend on the 
diffusion rates. 



II) A ^ A'. 

For this case one arrives at 



du = [(A — S) — (X' - 6)u]u 

da {\ + 6u)a ' 1 ' 



3 



which after integration gives, 

; 2 _ 5.' \"i /rx_\'\ 

-5+\ 

'u ' \a ' 



1 _ (X'-S)u (^ 2 -VA)/(5-A') 
1 A-<5 



(\'—S)u 



\-8 



(-)- A = - ) (io) 



Now, it is easy to obtain large time behavior of the average densities. Generally, 
there are three cases, 

III) 6> A, A' 

Depending on the initial average densities, the large time behavior of the average 
densities ratio is u(t) — > (b ~ a s / x ) or u — > oo (6 - a A '/ 5 ), which means that 
one kind of species decays faster. 

11.2) S< A, A' 

, [(i-A)(A'-i)] 

Defining e := w — -p-^ , at large times e(t) ~ a i**'-^] . In this case two kind 
of species decays with the same rate. 

11.3) A < 6 < X' 

At large times the average densities ratio u(t) — > 0, and 6 ~ a 5 / A . 

As it is seen, in the MF approximation only for a special choice of parameters, 
which is independent of diffusion rates, the two types of species decay with the 
same rate. The case with A = A' < 5, and Da = Db has been considered 
in • Using field theoretic renormalization group analysis, it is shown that 
in one-dimension both type of species decay with the same rate. Monte Carlo 
data also supports the field theory predictions in the one-dimensional model. 

3 The Cluster Mean Field Approximation 

Now, we want to use cluster mean-field approximation. If the diffusion rate 
for both type of species is the same, N = 2 cluster mean-field approximation 
gives the same value for the decay rates for both type of species, even if there 
is an imbalance in the initial average densities. If two type of species diffuse 
with different rates, irrespective of the initial values, at large times particles with 
greater diffusion rates decay more rapidly. For the nearest-neighbor interactions, 
the evolution equation of fc-point functions {nifio, ■ ■ ■ nu) contains at most (fc+1)- 
point functions. So, generally this set of equations will be a hierarchy, which 
can not be solved exactly. One way to overcome this difficulty is to impose 
constraints on the reaction rates that leads to disappearance of (k + l)-point 
functions from the evolution equation of fc-point functions. This method has 
been used to calculate some correlators exactly in 1 . Another possible way is 
to use the cluster approximation. In the fc-site cluster approximation, the set of 
evolution equations truncates and one encounters with a closed set of equations 
which may be solvable, at least numerically. Any how, for a two-site cluster 
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approximation, a three site joint probability for a sequence of nearest-neighbor 
sites is approximated by 

P(A, B, C) = P(A | B, C)P(B, C) c P[AB ^^ C) . (11) 

where P(A \ B) is the conditional probability. In the mean-field approximation 
there are three variables, (A), (B), and (0), only two of them are generally 
independent. In the two-site cluster approximation, or pair approximation, the 
variables are (A), (B), (0), (A0), (50), (00), among them there are six 
independent variables which we choose to be (A), (B), (AA),(BB), (AB), and 
(BA). In fact in the pair approximation, besides the average densities the two- 
point functions can also be obtained. The equation of motion for the average 
densities when Da = D b =: D are 

= -MBB) - t{AB) - 5 -(BA). (12) 

Similar to the mean-field approximation, the diffusion rates do not appear in the 
evolution equations of the average densities. But in fact they affect the average 
densities through the evolution equations of two-point functions, which are, 

^ = - \{AA) MAAA) S -(AAB) "-(BAA) 

- D(AA%) - D($AA) + 2D(A®A), (13) 

« = - \{BB) - X(BBB) - i(BBA) - '-(ABB) 

- D(BB$) - D($BB) + 2D(B$B), (14) 

^1 = S -(AB) S -(BAB) S -(ABA) D(AB%) 

- D(®AB) +2D(A®B), (15) 

= - \{BA) S -(BAB) S -(ABA) D(BA%) 

- D(®BA) + 2D(B$A). (16) 

To solve these equations in the cluster approximation, one should first approx- 
imate three-point functions and then all the equations should be expressed in 
terms of independent variables. For example (AB$) can be written as 

(^>,™ (17, 

and then using probability conservation (B%) should be expanded, 

(50) - (B) - (BA) - (BB). (18) 
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3.1 D A = D B 



Figure 1 and 2 show results for (A), (B), and the density ratios u(t) := (B) / (A) 
obtained using numerical solutions of equations 1(111150 . As it is seen both 
types of species decay with the same rate irrespective of equality or inequality 
of reaction rates A and A'. In the MF approach, K{t) := (BA)/(AA)is not 
an independent quantity and is (B)/(A). But in the CMF approach it is an 
independent one and the numerical result obtained for it is plotted in figure 3. 
As it is seen, in the CMF approximation it approaches a constant value at large 
times, means that both (AA) and (BA) decay with the same rate. Equality of 
their decay rates is independent of equality or inequality of reaction rates A and 
A'. 



3.2 D A D B 

As MF equations are independent of diffusion rates, their solutions remain un- 
altered. But in the pair approximation, only equations 1(12(1 remains unaltered. 
The diffusion rate D in 1(13114(1 should be changed properly to D A or D B , and 
the equations 1(15116(1 become. 



^1 = 6 -(AB) S -(BAB) S -(ABA) D B (AB<b) 

-D A {$AB) + {D A + D B ){A$B), (19) 



ffi^ = - \{BA) 5 -{BAB) - 5 -(ABA) - D A (BA9) 

-D B (d)BA} + (D A + D B )(B9A). (20) 

These set of equations has been solved numerically, and the numerical results for 
the average densities has been plotted in figures 4 and 5. As it is seen, at large 
times, species with greater diffusion rate dies out faster. If species with greater 
diffusion rate are majority initially, there is a cross over, as it is seen from figure 
4. The general behavior of the average densities ratio is independent of A, A', 
and 5, and the general form of large time behavior is determined by the diffusion 
rates. See figure 6. The numerical results for K(t) have been summarized in 
figure 7. 
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Figure captions 



figure 1- Average densities (A), and (B) as a function of time. The rates are 
D A = D B = 1, A = A' = 1000, and 5 = 3000. 

figure 2- Ratio of average densities, u(t) = (A) /(B), as a function of time. For 
the dashed line, the rates are D A = D B = 1, A = 1000, A' = 500, and 6 = 3000. 

figure 3- K(t) = (BA) / (AA) , as a function of time. 

figure 4- Average densities as a function of time. The rates are D A = 0.1, 
£> B = 1. 

figure 5- Average densities as a function of time. The rates are D A = 1, Db = 
0.1. 

figure 6- u(t) = (B)/(A), as a function of time. The rates are A = 2400, 
A' = 1000, 5 = 2500, and the greater diffusion rate is 1, and the smaller one, is 
0.1. 

figure 7- K(t) = (BA)/(AA), as a function of time. The rates are A = 250, 
A' = 100, 5 = 300, and the greater diffusion rate is 1, and the smaller one, is 
0.1. 
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Figure 1 
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